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We present the results of our analysis of the frequencies of galaxies with tidal tails and M51-type galaxies 
in several deep fields of the Hubble Space Telescope (HDF-N, HDF-S, HUDF, GOODS, GEMS). In total, 
we have found about seven hundred interacting galaxies at redshifts z < 1.5 in these fields. At 2 < 0.7, 
the observed space densities of galaxies with tidal structures and M51-type galaxies have been found to 
increase as cx (1 + z)™, where m « 2.6. According to our estimates, over the last 6-7 Gyr, i.e., at z < 0.7, 
about a third of the galaxies with M{B) < —18™ must have undergone strong gravitational perturbations 
and mergers and ~ 1/10-1/5 of the galaxies have accreted relatively low-mass nearby satellites typical 
of M51-type galaxies. The possible decrease in the timescale on which a distant galaxy appears peculiar 
with growing z can increase considerably the estimated rate of mergers. 



C^l. Introduction 
O 

)H The idea about the growth of stellar systems inside 
^ hierarchically merging dark haloes provides a basis 
I— ifor the present-day models of galaxy formation and 
I evolution (see, e.g., White and Rees 1978; for a re- 
J> view, see the book by Mo et al. 2010). The process 
0^ of halo mergers can be well described in terms of 
piodel calculations (Fakhouri and Ma 2008 and refer- 
ences therein); however, the connection of this process 
with the mergers of actually observed galaxies is not 
yet understood well enough (see, e.g., Kitzbichler and 
^ White 2008). 

^1 As a rule, two approaches are used when the galaxy 
merger rate is studied observationally. The first is 
5_i based on analysis of the frequency of close galaxy 
pairs at various z (Zepf and Koo 1989; Kartaltepe et 
al. 2007; de Ravel et al. 2011; and references therein). 
The weak points of this method are the small num- 
ber of spectroscopically confirmed pairs at high z and 
the assumption that all these pairs merge on a certain 
timescale (Kitzbichler and White 2008). 

The second (morphological) approach is based 
on the frequency statistics of signatures of re- 
cent interactions and mergers in galaxies: struc- 
ture distortions (Conselice et al. 2003 and refer- 
ences therein) , the presence of tidal tails and bridges 
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(Reshetnikov 2000a, 2000b; Bridge et al. 2010), polar 
rings (Reshetnikov 1997; Reshetnikov and Dettmar 
2007), coUisional rings (Lavery et al. 2004), etc. The 
main problem of this method is that the basic trac- 
ers of dynamical perturbations of galaxies, as a rule, 
have a low surface brightness and they are difficult to 
observe at high z. However, as was shown by Hibbard 
and Vacca (1997), such structures must be visible at 
least up to 2; ~ 1 on deep images of the Hubble Space 
Telescope (HST). The search for and statistics of tidal 
structures in galaxies in the HDF-N and HDF-S deep 
fields confirmed this conclusion (Reshetnikov 2000a, 
2000b). Furthermore, as in the case of studying merg- 
ing pairs, the connection of the dark halo merger rate 
derived from models with the complex processes of 
galaxy interactions and mergers acting on different 
mass and timescales is ambiguous. On the other hand, 
the advantage of the morphological approach is that 
very large samples of galaxies can be produced and 
investigated using it. 

The results of the two approaches do not yet agree 
quite well quantitatively; however, the general con- 
clusion is beyond question: the fraction of interact- 
ing and merging galaxies actually grows with z, as 
follows from theoretical expectations. At z < 1.5, 
this growth is commonly described by a power law, 
(l + z)"^, where the value of m from the data of differ- 
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ent authors is, in most cases, within the range from 
2 to 4 (see, e.g., Table 2 in Kartaltepe et al. 2007). 

In this paper, we analyze the frequency statistics of 
galaxies with tidal structures and M51-type galaxies 
in several Hubble deep fields. All numerical values in 
the paper are given for the cosmological model with a 
Hubble constant of 70 km s~^ Mpc~^ and = 0.3, 
Qa = 0.7. 

2. Galaxies with tidal structures at z ^ 0.7 
2.1. The sample of galaxies 

The sample is based on the catalog of interacting 
galaxies in several Hubble deep fields (Mohamed and 
Reshetnikov 2011). This catalog was compiled on 
the basis of a visual classification of galaxies in the 
F814W (HDF-N and HDF-S fields), F775W (HUDF), 
and F850LP (GOODS and GEMS) filters. At z ~ 1, 
these filters roughly correspond to the B band in the 
reference frame associated with the galaxies them- 
selves. Our sample includes a total of 689 galaxies 
with tidal tails (the bridges have, on average, a lower 
surface brightness and they were not considered in 
our work) at redshifts z < 1.5. 

Figure 1 shows the z distribution of sample galax- 
ies and their positions on the "rest-frame absolute 
magnitude - redshift" diagram. The dashed lines 
in the upper panel show the fits of the observed 
distributions by the empirical formula dN/dz oc 

exp( — (^r/zc)^'^ proposed for describing the obser- 
vational data for magnitude-limited samples with me- 
dian redshift z^ = 1.4:12zc (Baugh and Efstathiou 
1993). As we see from the figure, the actual distribu- 
tions agrees satisfactorily with the expected ones. 

Figure 1 (bottom part) illustrates the obvious ob- 
servational selection effect: we predominantly choose 
the brightest galaxies among more distant objects. 
The standard way of avoiding this selection is to con- 
sider only bright galaxies at various z. Below, we will 
study the statistics of interacting galaxies for two 
subsamples: (1) for objects with absolute B magni- 
tude from -18™ to -20™" and (2) for all galaxies with 
M{B) <-18™ (the stellar mass limit M> 4 • 10^ Mq 
roughly corresponds to this luminosity restriction). 

To estimate the completeness of our sample, we 
considered the differential counts of galaxies in 
comparison with those in the VIRMOS VLT F02 
deep field (McCracken et al. 2003). According to 
McCracken et al. (2003), the differential counts for 
galaxies brighter than 24"^ in the / band have a slope 
of 0.34ib0.02. (This result agrees well with the data 
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Fig. 1. Top: Redshift distributions of the samples of 
galaxies with tidal tails (red histogram) and M51- 
type objects (blue histogram). The dashed lines in- 
dicate the expected distributions for Zm = 0.64. 
Bottom: The absolute magnitude - redshift relation 
for galaxies with tidal tails (red circles) and M51-type 
galaxies (blue diamonds). 



on other deep fields as well; e.g., Metcalfe et al. 2001). 
The objects in our sample of interacting galaxies fol- 
low this slope up to / ~ 21.5 and then begin to devi- 
ate greatly from it. The median redshift for galaxies 
with / ~ 21.5 is z — 0.67 (this value is also close to 
that found from the analytical fit shown in Fig. 1) and 
below we take precisely this redshift complete- 
ness limit for the sample of interacting galaxies (both 
for galaxies with tidal tails and M51-type objects). 

2.2. Estimating the evolution rate 

To study the evolution of the frequency of galaxies 
with tidal structures, we used the same approach as 
that in Reshetnikov (2000a, 2000b). The meaning of 
this method is that, having fixed the space density of 
objects of some type at z = 0, the expected number 
of such galaxies is estimated within the selected field 
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in a given z range for different space density evolution 
laws. By comparing the actual and expected numbers 
of objects, we can obtain a constraint on the density 
evolution law. Below, we assume that the space den- 
sity of galaxies varies with z as n{z) = no x (1 + z)™ 
and estimate the exponent m. 

An important stage in estimating the evolution 
rate is to determine the local density of galaxies ng. 
Unfortunately, for galaxies with tidal structures, this 
quantity is known relatively poorly. In addition, it de- 
pends on the brightness level at which the structure 
is distinguished: for example, at a brightness level 
of ~ 28'"/n", more than 10% of spiral galaxies ex- 
hibit various kinds of external structures (Miskolczi 
et al. 2011). In our paper, we used the results of Nair 
and Abraham (2010), who performed a detailed vi- 
sual classification of approximately 14000 galaxies 
from the Sloan Digital Sky Survey; they also provided 
the frequency statistics of galaxies with tidal tails 
(see Table 4 in their paper). According to Nair and 
Abraham, the visually distinguishable tidal tails are 
seen in 301 of 14034 galaxies, i.e., in about 2% of the 
galaxies. As the luminosity function of nearby galax- 
ies, we take the results of the 2dF survey (Norbcrg 
et al. 2002) and assume that the fraction of galaxies 
with tidal tails in any luminosity range is 0.02. 

When nearby and distant galaxies are compared, 
the possible evolution of their luminosity with z 
should be taken into account. Studies of the Tully- 
Fisher relation and the luminosity function of distant 
galaxies show that the galaxies at z ~ 1 were approx- 
imately 1™ brighter (see, e.g., Gabasch et al. 2004; 
Miller et al. 2011). Therefore, to properly compare 
the space densities of galaxies at different 2, this ef- 
fect should be taken into account. As a model for the 
luminosity evolution of spiral galaxies, we used the 
results of Bicker et al. (2004). In the redshift range 
of interest [z < 0.7), the predictions of the model 
by Bicker et al. for the luminosity evolution of Sb-Sc 
spirals (AM(B) oc (1"^ - l."^2) x z) are close to the 
actually observed changes. 

Table 1 presents the final results of determining the 
exponent m from the statistics of galaxies with tidal 
tails within the deep fields wc considered (their to- 
tal angular area is = 6.83 x 10^^ sr) at z < 0.67. 
As the errors in m, the table provides the range of 
variation in the exponent in the case of a Poissonian 
error in the number of objects in the field (it\/iV); 
the total number of galaxies in a given luminosity 
range is given in parentheses ( A^) ; the first and second 
rows give the value of m in the absence of luminosity 



Table 1. Results of the determination of the param- 
eter m 



Sample 


-18" >Mb> -20" 


Mb < -18" 


Galaxies 
with tidal tails 

M51-type 
galaxies 


2.561°;^^ (198) 
2.nt°o\t (211) 

2.97t°i3 (27) 
2mlL°oil (24) 


2.93V^:lt (278) 
2.60t°:l^ (243) 

s-nt'oil (34) 

2.48t°-il (26) 



evolution (only the ^-correction was applied to the 
absolute magnitudes) and with its allowance. As we 
see from the table, allowance for the luminosity evo- 
lution changes the m estimates, but not too much, 
approximately by 10%. 

3. The frequency of M51-type galaxies 

M51-type galaxies are binary systems that consist, 
as a rule, of a bright spiral galaxy with a relatively 
low-mass satellite located near the end of one of its 
spiral arms (Klimanov and Reshetnikov 2001). In 
the nearby Universe, such systems are relatively rare 
(only about 0.3% of all galaxies can be attributed to 
this type; see Klimanov 2003). However, they are con- 
venient objects for studying a number of questions in 
the physics of galaxies (e.g., the spiral arms formed 
by tidal perturbations, the effects of satellites on the 
structure and star formation in galaxies). 

As in the case of galaxies with tidal tails, the sam- 
ple of M51-type galaxies is based on the catalog by 
Mohamed and Reshetnikov (2011); see the examples 
in Fig. 2. The sample includes 74 such binary sys- 
tems (the distributions of galaxies in z and M{B) 
are shown in Fig. 1). Figure 3 compares the distri- 
butions of nearby (32 binary systems; see Klimanov 
and Reshetnikov 2001) and distant galaxies in ratio 
of the observed luminosity of the satellite to the lu- 
minosity of the main galaxy [Lg/Lm) and in relative 
distance to the satellite {R/r). As we see from the 
figure, the distributions are, on the whole, similar. 
The corresponding means are {Ls/Lm) = 0.19ib0.21, 
{R/r) = 1.34±0.48 (nearby galaxies) and 0.22±0.18, 
1.26 lb 0.53 (distant galaxies). 

The local space density no of M51-type galaxies 
was found from their luminosity function. According 
to Klimanov (2003), the i?-band luminosity function 
can be described by a Schechter function with the 
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Fig. 2. Examples of M51-type galaxies in the GEMS field. The images were extracted from the frames in 
the F850LP filter. The image sizes for the first four galaxies are 5."4 x 5."4, 3."6 x 3."6 for 37254, 18" x 18" 
for 41271. The galaxy number and redshift according to the catalog of Wolf et al. (2004) are indicated in the 
lower left corner of each figure. 



following parameters: = 1.18 x 10 ^ Mpc ^, a = 
-1.30, M^B) = -20.45. Then, noi-lS"" > M{B) > 
-20"^] = 2.51 x 10-^ Mpc-3 and no[M{B) < 
-18"*] = 2.97 x 10-^ Mpc-'l The final estimates of 
the parameter m for two luminosity ranges are listed 
in Table 1. As in the case of galaxies with tidal tails, 
the first and second rows give the estimates without 
and with allowance for the luminosity evolution, re- 
spectively. 

4. Discussion 

4-1- Galaxies with tidal tails 

Let us consider how our data agree with the results 
of other works. Methodically, the closest work is the 
paper by Reshetnikov (2000b), in which it was found 
from the statistics of 25 galaxies with tidal tails in 
HDF-N and HDF-S that m = S-Gto^. To properly 
compare the results, we should take into account the 
fact that Reshetnikov (2000b) assumed the galaxies 
with tidal tails to account for 1% of the field galaxies. 
If this estimate is doubled (see Sect. 2.2), then the m 
estimate will decrease approximately by a factor of 
1.5 and will be consistent with our results. 



The results of several recent works, in which the 
morphological approach, i.e., the statistics of signa- 
tures of gravitational perturbations and mergers, was 
used to estimate m at z ~ 1 are summarized in 
Table 2 (see also references in these papers). The first, 
second, third, fourth, fifth, and sixth columns of this 
table give, respectively, the references, the type of 
objects ('m' for mergers, 'rg' for ring galaxies, 'ts' 
for tidal structures, 'tt' for tidal tails), the number 
of galaxies in the sample of interacting galaxies (if 
it is given in the original paper), the redshift range, 
the range of M{B) or another similar characteristic 
(the infrared luminosity, the mass of stars), and m. 
Our results on M51-type galaxies were not included 
in Table 2, because they should be compared with the 
statistics of close pairs. 

As we see from Table 2, the works on the statis- 
tics of tidal structures (this paper and the paper by 
Bridge et al. 2010) yield very similar results. Bridge 
et al. performed a visual classification of about 27 000 
galaxies with i < 22."* 2 in a sky region with an area 
of 2n° (the CFHTLS-Deep project). As a resuh, they 
identified more than a thousand interacting galaxies 
with z from 0.1 to 1.2. Comparing the number of 
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Reference 


Type of objects 


N 


Az 


AM{B) 


m 


Conselice et al. (2003) 


m 




< 1 


< -18'" 


2.5±0.3 


Lavery et al. (2004) 


rg 


25 


0.1-1.0 




5.2±0.7 


Bridge et al. (2007) 


m 


64 


0.2-1.3 


LiB,>5- 101°Lq 


2.12±0.93 


Kampczyk et al. (2007) 


m 


68 


< 0.7 


< -19" 


3.8±1.2 


Lotz et al. (2008) 


m 


402 


0.2-1.0 


< -18.94- 1.32 


0.23±1.03 


Conselice et al. (2009) 


m 




0.2-1.2 


> 101"Mq 


2.3±0.4 


Lopez-Sanjuan et al. (2009a) 


m 


25 


0.35-0.85 


< -20" 


2.9±0.8 


Lopez-Sanjuan et al. (2009b) 


m 


61 


0.2-1.1 


< -20" 


1.8±0.5 


Bridge et al. (2010) 


ts 


1075 


0.25-1.0 


> 10^-5 M0 


2.56±0.24 


This paper 


tt 


243 


< 0.67 


< -18" 






tt 


211 


< 0.67 


[-18", -20"] 


9 71 0.16 
^- ' -"^-0.17 
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Fig. 3. Top: Ratio of the luminosity of the satellite to 
the luminosity of the main galaxy for the samples of 
nearby (black solid line) and distant (red dashed line) 
M51-type galaxies. Bottom: Ratio of the distance to 
the satellite (R) to the radius of the main galaxy (r) 
for M51-type galaxies (the notation is the same). 



galaxies with tidal structures with the total number 
of objects in this z range, they constructed the red- 



shift dependence of the fraction of interacting galax- 
ies. The derived dependence is described by the ex- 
ponent m = 2.25 lb 0.24; however, if the least reliable 
intervals near the beginning and the end of the inves- 
tigated z range are eliminated, then m = 2.56 it 0.24 
(this value is given in Table 2). Note also that Bridge 
et al. (2010) considered galaxies with stellar masses 
> 10^'^ Mq (Table 2). Our approach to estimating 
m and the observational data used were completely 
different, but the results turned out to be close. 

The paper by Lavery et al. (2004) on collisional 
ring galaxies gives a very large value of m. However, 
the sample of galaxies they used is small. In addi- 
tion, the type of interaction itself in which one of the 
galaxies passes through the center of the other at a 
high velocity is so rare that the local space density 
of such objects and, consequently, the value of m are 
known poorly. 

The statistics of merging galaxies (mergers) also, 
on the whole, is indicative of relatively large (w 2 — 4) 
m (Table 2). The only exception is the paper by 
Lotz et al. (2008); however, there are several fac- 
tors (e.g., a strong dependence of the result on the 
data in two bins with the lowest z) that may ex- 
plain this discrepancy (for a discussion, see Bridge 
et al. 2010). As a result, we may conclude that the 
works of recent years give a relatively consistent pic- 
ture of the growth of the frequency of morphological 
signatures of interactions and mergers up to z ~ 1. 
Most of the results agree, within the error limits, with 
m ~ 2.6. The statistics of binary systems, including 
M51-type galaxies (this paper), also yields similar re- 
sults, m ~ 3 (Le Fevre et al. 2000; Kartaltepe et al. 
2007; de Ravel et al. 2011). 
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The conversion from the observed fraction of inter- 
acting and merging galaxies to the galaxy merger rate 
{Rmg) is usually accomplished through the timescale 
tv on which the galaxies appear morphologically pe- 
culiar (e.g., exhibit distinguishable tidal structures): 
Rmg{z) = 5{z)/ty, where S{z) is the fraction of inter- 
acting galaxies at redsliift z. For galaxies with tidal 
tails, we assumed that Sq = 0.02; at m = 2.6 for 
galaxies with M{B) < -18™ (Table 2), we then ob- 
tain S{z = 0.7) = 0.08. As the time during which the 
galaxies exhibit noticeable tidal tails, we take, accord- 
ing to Bridge et al. (2010), ty— 0.8 Gyr. (Of course, 
this estimate is very uncertain, because it depends, in 
particular, on the mass ratio of the interacting galax- 
ies, the gas fraction in them, the method of identify- 
ing tidal structures, etc.) For the galaxy merger rate 
at z = 0.7, we then obtain Rmg ~ 0.1 mergers per 
galaxy during 1 Gyr. This estimate agrees well with 
the data of other authors obtained by different meth- 
ods (sec Fig. 11 in Bridge et al. 2010). If we pass to 
a unit volume, then the merger rate is ~ 10~^/ (Gyr 
X Mpc^). 

Using the above numbers, we can estimate the av- 
eraged history of interactions for the galaxies of our 
sample. Integrating the fraction of interacting galax- 
ies at different z by taking into account the timescale 
ty (see Eq. (3) in Bridge et al. 2010), we found that 
a typical galaxy with M{B) < —18™ in the z range 
from 0.7 to 0.0 (a time of 6-7 Gyr corresponds to this 
range) underwent ~ 0.35 mergers or close encounters 
accompanied by the formation of extended tidal tails. 
In other words, a third of bright galaxies have under- 
gone strong gravitational perturbations and mergers 
in the last 6-7 Gyr. 

In the above simple estimates, there is one, previ- 
ously almost undiscussed uncertainty: the timescale 
ty can depend on z. As was shown by Mohamed et 
al. (2011), the tidal structures in distant galaxies, on 
average, appear shorter than those in nearby objects. 
One of the possible causes is observational selection 
due to the cosmological brightness dimming and the 
influence of the fc-correction. As a result of this selec- 
tion, in distant galaxies we predominantly observe a 
relatively early evolutionary stage of the tails, when 
they have a high surface brightness, while in nearby 
galaxies we see, on average, "older" and longer struc- 
tures. 

To quantitatively estimate the influence of this ef- 
fect on Rmg, we used the results of the calculations 
by Mihos (1995). Mihos showed that the tidal tails 
should be visible in the HST exposures only for ~ 150 



Myr at z = 1 and for about 350 Myr at z = 0.4. 
Assuming that ty = 800 Myr for z = (see above), 
we can obtain the following simple approximation for 
the z dependence of ty: ty = aexp{—z/b) + c, where 
a = 0.7, b = 0.4, c = 0.095, and the timescale ty is 
in Gyr. For z = 0.7, it follows from this approxima- 
tion that ty = 0.22 Gyr and, accordingly, Rmg ~ 0.36 
mergers per galaxy during 1 Gyr. For the full history 
of mergers between z = 0.7 and z = 0.0, we obtain 
an estimate of ~0.75, i.e., 3/4 of all galaxies with 
M{B) < —18"* must have undergone mergers in the 
last 6-7 Gyr. 

4. -2. M51-type galaxies 

Let us now consider M51-typc galaxies. As we showed 
in this paper, their space density at z < 0.7 evolves 
with m = 2.5 — 2.7 (Table 1). Assuming that m = 2.6 
and 6(j — 0.003 (Klimanov 2003), we can estimate 
that the relative fraction of such galaxies is 0.012 
at z = 0.7 and 0.014 at z = 0.8. Recently, Lopez- 
Sanjuan et al. (2011) presented the results of direct 
counts of faint satellites (with luminosities from 0.1 
to 0.25 of the luminosity of the central galaxy) for 
various distances from the main galaxy for galax- 
ies with M{B) < -20"* at z = 0.5 and z = 0.8 
based on the VIMOS VLT data (VVDS-Deep sur- 
vey). Extrapolating these counts to 15 kpc (most of 
the satellites in our sample of M51-type galaxies are 
within this distance), we can estimate that the frac- 
tion of such galaxies at 2; = 0.8 is «0.02. Given the 
large errors in such estimates (different methods of 
identifying objects, different constraints on their lu- 
minosity, etc.), the agreement between the estimates 
of the fraction of galaxies with relatively low-mass 
satellites at z = 0.8 may be recognized to be satisfac- 
tory. 

The satellites near the outer boundary of the stellar 
disks in the main galaxies must be rapidly swallowed. 
The characteristic merger timescale is ~ 0.2 — 0.4 Gyr 
(see, e.g., Lotz et al., 2010). Consequently, the merger 
rate estimate for M51-type galaxies at 2: = 0.7 is 
Rmg Si 0.03 - 0.06 per galaxy brighter than M{B) = 
—18™ in 1 Gyr. Lopez-Sanjuan et al. (2011) give a 
similar value: the minor merger rate at z = 0.5 — 0.8 
is 0.0045-0.034 (Table 8 in their paper). The total ex- 
pected number of mergers per galaxy at 2; < 0.7 is 
0.1-0.2. Therefore, M51-type galaxies make a notice- 
able and non-negligible (compared to major mergers) 
contribution to the galaxy merger rate at 2 = 0.7. Of 
course, this conclusion is very sensitive to the satel- 
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lite merger timescale. If, for instance, it is close to 1 
Gyr (Kitzbichler and White 2008), then the merger 
rate through M51-type galaxies decreases by several 
times. 

5. Conclusions 

Based on a large sample of distant galaxies with tidal 
tails and M51-type galaxies (spiral galaxies with a 
relatively low-mass satellite located near the end of 
one of the spiral arms), we estimated the evolution 
of the space densities of objects of these types up to 
z = 0.7. It turned out that their observed densities 
increases with redshift approximately as (1 -|- z)^'^. 

At z = 0.7, the merger rate leading to the for- 
mation of extended tidal tails is ~ 0.1 per galaxy 
brighter than M{B) = —18™" in 1 Gyr. The corre- 
sponding merger rate for M51-type galaxies is approx- 
imately a factor of 2-3 lower. 

In the last 6-7 Gyr, i.e., at 2; < 0.7, about a third 
of the galaxies with M{B) < —18™ must have under- 
gone strong gravitational perturbations and mergers; 
~ 1/10-1/5 of the galaxies swallowed nearby satellites 
with Ls/Lm 0.1 — 0.2. Such processes are capable 
of radically changing the characteristics of galaxies 
(Toomre, 1977) and stimulating the processes of star 
formation and nonthermal activity of the nuclei (see, 
e.g.. Keel et al. 1985). 

The estimates of the galaxy merger rate depend 
strongly on the adopted timescale on which they ap- 
pear peculiar (i^). For instance, allowance for the pos- 
sible redshift dependence of (the identification time 
of tidal structures decreases with increasing z) can in- 
crease the above merger rates by several times. 

The high galaxy merger rates found from obser- 
vations of distant objects clearly suggest that grav- 
itational interactions and mergers were among the 
most important processes that determined the in- 
dividual properties of the galaxies surrounding us. 
Furthermore, they affect the evolution of the "lumi- 
nous" matter in the Universe whole, because 
mergers change the luminosity function of galaxies, 
the luminosity density produced by them, and other 
characteristics. 
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